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This paper deals with non-Markovian behaviour in atomic systems coupled to a structured 
reservoir of quantum EM field modes, with particular relevance to atoms interacting with the field 
in high Q cavities or photonic band gap materials. In cases such as the former, we show that the 
pseudo mode theory for single quantum reservoir excitations can be obtained by applying the Fano 
diagonalisation method to a system in which the atomic transitions are coupled to a discrete set of 
(cavity) quasimodes, which in turn are coupled to a continuum set of (external) quasimodes with 
slowly varying coupling constants and continuum mode density. Each pseudomode can be identified 
with a discrete quasimode, which gives structure to the actual reservoir of true modes via the 
expressions for the equivalent atom-true mode coupling constants. The quasimode theory enables 
cases of multiple excitation of the reservoir to now be treated via Markovian master equations for 
the atom-discrete quasimode system. Applications of the theory to one, two and many discrete 
quasimodes are made. For a simple photonic band gap model, where the reservoir structure is 
associated with the true mode density rather than the coupling constants, the single quantum 
excitation case appears to be equivalent to a case with two discrete quasimodes. 

PACS numbers; 42.50.-p,42.70.Qs,42.50.Lc 



I. INTRODUCTION 

The quantum behaviour of a small system coupled to a large one has been the subject of many studies since 
quantum theory was first formulated. The small system is usually of primary interest and generally microscopic 
(atom, nucleus, molecule or small collections of these) but currently systems of a more macroscopic nature (Bose 
condensate, superconductor, quantum computer) are being studied. The large system is invariably macroscopic in 
nature (free space or universe modes of the EM field, lattice modes in a solid, collider atoms in a gas) and is of less 
interest in its own right, being primarily of relevance as a reservoir or bath affecting the small system in terms of 
relaxation and noise processes. The large system is often a model for the entire external environment surrounding the 
small system. Changes in the small system states (described in terms of its density operator) can be divided into two 
sorts — effects on the state populations (energy loss or gain) or effects on the state coherences (decoherence or induced 
coherence). Equivalently, quantum information (described via the von Neumann entropy) would be lost or gained 
due to the interaction with the environment, and its loss is generally associated with decoherence. Interestingly, as 
the small system becomes larger or occupies states that are more classical the time scale for decoherence can become 
much smaller than that for energy loss. This is of special interest in quantum information processing [|l]-^ where 
the small system is a collection of qbits making up a quantum computer weakly coupled to the outside world, or in 
measurement theory ^-Q], where the small system is a micro system being measured coupled to an apparatus (or 
pointer) that registers the results. For quantum computers it is desirable that decoherence is negligible during the 
overall computation time |^ (otherwise error correction methods have to be incorporated, and this is costly in terms 
of processing time), whereas in measurement theory environment induced decoherence is responsible for the 
density operator becoming diagonal in the pointer basis (otherwise a macroscopic superposition of pointer readings 
would result). 

A standard method for describing the reservoir effects on the small system is based on the Born-Markoff master 
equation for the system density operator [p^|-p^ . This depends on the correlation time for the reservoir (as determined 
from the behaviour of two time correlation functions for pairs of reservoir operators involved in the system reservoir 
interaction) being very short compared to that of the relaxation and noise processes of the system. In general 
terms, the more slowly varying the coupling constants for this interaction or the density of reservoir states are with 
reservoir frequencies, the shorter the correlation time will be. For many situations in the fields of quantum optics, 
NMR, solid state physics the Born-Markoff master equation provided an accurate description of the physics for the 
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system. Elaborations or variants of the method such as quantum state trajectories 15|, Fokker-Planck or c- 
number Langevin equations quantum Langevin equations are also used. Sometimes an apparently 

non-Markovian problem can be converted to a Markovian one by a more suitable treatment of the internal system 
interactions (for example, the use of dressed atom states [ ^t|Ji8| for treating driven atoms in narrow band squeezed 
vacuum fields 

However, situations are now being studied in which the standard Born-Markoff approach is no longer appropriate, 
since the reservoir correlation times are too long for the time scales of interest. A structured rather than flat reservoir 
situation applies ||2l] , ^2[| . This includes cases where the reservoir coupling constants vary significantly with frequency, 
such as the interaction of atom(s) or quantum dots with light in high Q cavities j23|, including microcavities (see, 
for example, [p4|-p8|) and microspheres (see, for example, |29 30|]) . Cases such as an atom (or many atoms — super- 
radiance) interacting with light in photonic band gap materials [ pl| , ^ ^ 33 1 , [ p4|j35| , where the reservoir mode densities 
that have gaps and non-analytic behaviour near the band gap edges, also occur. Also, quantum feedback situations 
p6| , p7| can involve significant time delays in the feedback circuit, and thus result in non-Markovian dynamics for the 
system itself. Furthermore, systems with several degrees of freedom, such as in quantum measurements (for example, 
the Stcrn-Gerlach experiment) could involve situations where the decoherence times associated with some degrees of 
freedom (such as the position of the atomic spin) could become so short that the Markoff condition might no longer 
be valid, and the effects of such non-Markovian relaxation on the decoherence times associated with more important 
degrees of freedom (atom spin states) would be of interest. 

A number of methods for treating non-Markovian processes have been developed. Apart from direct numerical 
simulations [Q, these include the Zwanzig-Nakajima non-Markovian master equation and its extensions 
the time-convolutionless projection operator master equation Ea], Heisenberg equations of motion p3|,p3[, stochastic 



wave function methods for non-Markovian processes |44-4£|, methods based on the essential states approximation 
or resolvent operators [^,^ 50 1, the pseudo mode approach |5^j5^, Fano diagonalisation p3| , p^j54| and the sudden 
decoherence approximation~[55[ . Of these methods the last four are easier to apply and give more physical insight into 
what is happening. However the essential states method is difficult to apply in all but the simplest situations, since 
the set of coupled amplitude equations becomes unwieldy and it is difficult to solve — as in the case where multiple 
excitations of the reservoir are involved. The pseudo mode approach is based on the idea of enlarging the system to 
include part of the reservoir (the pseudo mode — which could be bosonic or fermionic depending on the case) thereby 
forming a bigger system in which the Markoff approximation now applies when the coupling to the remainder of the 
reservoir is treated. At present the pseudomode method is also restricted to single reservoir excitation cases. The 
Fano diagonalisation method relates the causes of non-Markovian effects to various underlying features (such as the 
presence of bound states for treating atom lasers), and is closely related to the pseudo mode method. The sudden 
decoherence method enables decoherence effects on time scales short compared to system Bohr periods to be treated 
simply via ignoring the system Hamiltonian. 

This paper deals with the relationship between the current pseudo mode method for single quantum reservoir 
excitations and the Fano diagonalisation method for situations where the reservoir structure is due to the presence 
of a discrete, system of (quasi) modes which are coupled to other continuum (quasi) modes. This important case 
applies to atomic systems coupled to the quantum EM field in high Q resonant cavities, such as microspheres or 
microcavities. It is shown that the pseudo mode method for single quantum excitations of the structured reservoir 
can be obtained by applying the Fano diagonalisation method to a system featuring a set of discrete quasi modes 
p6| , ^ together with a set of continuum quasi modes, whose mode density is slowly varying. The structured reservoir 
of true modes (S^js^ is thus replaced by the quasi modes. The interaction between the discrete and continuum quasi 
modes is treated in the rotating- wave approximation and assuming slowly varying coupling constants ||60| , |6l|| . The 
atomic system is assumed to be only coupled to the discrete quasi modes. The density of continuum quasi modes is 
explicitly included in the model. Although the behavior of the atomic system itself is non-Markovian, the enlarged 
system obtained by combining the discrete quasi modes with the atomic system now exhibits Markovian dynamics. 
The discrete quasi modes are identified as pseudo modes. The continuum quasi modes are identified as the flat 
reservoir to which the enlarged Markovian system is coupled. Explicit expressions for the atom-true modes coupling 
constants are obtained, exhibiting the rapidly varying frequency dependence characteristic of structured reservoirs. 
At present the treatment is restricted to cases where threshold and band gap effects are unimportant, but may be 
applicable to two-dimensional photonic band gap materials. However, the problem of treating multiple excitation 
processes for certain types of structured reservoirs can now be treated via the quasimode theory, since the Markovian 
master equation for the atom-quasi mode system applies for cases involving multilevel atoms or cases of several excited 
two-level atoms. Further extensions of the treatment to allow for atomic systems driven by single mode external laser 
fields are also possible, with the original atomic system being replaced by the dressed atom. 

The plan of the paper is as follows. In Section || the key features of pseudo mode theory are outlined. Section 
Ip presents the Fano diagonalisation theory for the quasi mode system, with details covered in the Appendices 
^and |c[ In Section IV specific cases such as one or two discrete quasi modes or where the variation of discrete- 
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continuum coupling constants can be ignored are examined, 
and reservoir structure functions in these situations. Section 
plus discrete quasi modes system 



iving results for the atom-true mode coupling constants 
contains the Markovian master equation for the atom 



Section V] briefly examines the situation where coupling constants and mode 
densities are not slowly varying. Conclusions and comments are set out in Section VI] . 



II. PSEUDOMODE THEORY 

The simplest case to which pseudomode theory |^ can be applied is that of a two-level atom coupled to the modes 
of the quantum electromagnetic (EM) field — which constitutes the structured reservoir. Only one photon excitation 
processes will occur. However, the formalism would also apply to any spin 1/2 fermion system coupled to a bath of 
bosonic oscillators. The Hamiltonian is given in the rotating wave approximation by: 

H — fiijjxa\ax -I- 2 ^'^i {^^^^ ~ o" o'^) 

A 

+ E(^5A«A<T++h.C.) (1) 
A 

where &^,ct~ are the usual atomic spin operators, a\, a| are the annihilation, creation operators for the mode A of 
the field, oji is the atomic transition frequency, wa is the mode frequency and g\ are the coupling constants. This 
system is illustrated in Fig. |l|. 

To describe a one photon excitation process, the initial condition is the atom excited and no photons present in the 
field. Hence the initial Schrodinger picture state vector is: 

|*(0)) = |1)|...0,...) (2) 

In the essential states approach the state vector at later time t will be a superposition of the initial state and states 
with the atom in its lower state and one photon in various modes A. With ci, c\ defining the complex amplitudes for 
the atomic excited state and the one photon states, the state vector is: 

|*(t))=ci(Oe— ^*|l)|...0....)+^c.We— ^*|0)|...U...> (3) 

A 

Substitution into the time dependent Schrodinger equation leads to the following coupled complex amplitude equa- 
tions: 

■ d ~ * ~iAxt~ 

«^ci = ^5;,e ^ CA 

*^5a = (4) 
at 

where A a = uj\ — uji are detunings. 

Formally eliminating the amplitudes for the one photon states enables an integro-differential equation for the excited 
atomic amplitude to be derived. This is: 

^ci(i) = -^*drG(r)5i(t-r) (5) 

and involves a kernel G given by 

G(r)=5]|5Ape-^^^ 

A 

2„— iAvT 



dLuxp{cox)\gx\'e-''^^^ (6) 

The mode density p{lu\) is introduced after replacing the sum over A by an integral over the mode angular frequency 
ujx. 

It is apparent from Eq.(^ that the behaviour of the atomic system only depends on the reservoir structure function 
D{lli\) for this single quantum excitation case defined by 
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Px\gxf ^^Diiox) (7) 

where a transition strength CI is introduced to normahse D so that its integral gives 2tt. The transition strength CI is 
given by: 

= J duxpiLux)\gx\'' (8) 

The reservoir structure function D[ijJx) enables us to describe the various types of reservoir to which the atomic system 
is coupled. If D is slowly varying as a function of lox then the reservoir is 'flat', whilst 'structured' reservoirs are where 
D varies more rapidly, as seen in Figure There are of course two factors involved in determining the behaviour of 
D — the mode density p{ujx) and the coupling constant via \gx\'^- Either or both can determine how structured the 
reservoir is. Photonic band gap materials are characterised by mode densities which are actually zero over the gaps 
in the allowed mode frequencies, and which have non analytic behaviours near the edges of the band gaps. All mode 
densities are zero for negative lox, so threshold effects are possible. In cavity QED situations, such as for microspheres 
and other high Q cavities, the coupling constant varies significantly near the cavity resonant frequencies, so in these 
cases it is the coupling constant that gives structure to the reservoir. For the present it will be assumed that (apart 
from simple poles) D is analytic in the lower half complex lox plane, and any other non analytic features can be 
disregarded. It is recognized of course that this restriction places a limit on the range of applicability of the theory, 
though it may be possible to extend this range by representing the actual function D{ujx) in an approximate form 
that satisfies the analyticity requirements. In addition (in order to calculate the contour integrals) it will be assumed 
that D tends to zero at least as fast as 1/|wa| as \ujx\ tends to infinity. 

Based on the above assumption regarding the reservoir structure function D, the kernel G may be evaluated in 
terms of the poles and residues of D{lox) in the lower half complex lox plane. It is assumed that these simple poles 
can be enumerated. The poles are located at zi, Z2, ....z;... and their residues are ri,r2, ...r;.... The pole zi may be 
expressed in terms of a real angular frequency uji and a width factor F; via zi = loi — iTi/2. Contour integration 
methods show that the sum of the residues equals i. The kernel is obtained in the form: 

G(t) = ne-'^"^'''^^ (9) 

I 

The integro-differential equation (0) for the excited atomic amplitude involves a convolution integral on the right 



hand side and may be solved using Laplace transform methods. The atomic behaviour obtained is well known ^8 59 
and will not be rederived here. It is found that there are two regimes, depending on the ratio of the transition strength 
to typical width factors. These are: (a) a strong coupling regime with Non-Markovian atomic dynamics, and occuring 
when 17 3> F (b) a weak coupling regime with Markovian atomic dynamics, occuring when il <C F. 

The pseudomode approach continues by considering poles of the reservoir structure function D(ujx) in the lower half 
complex ujx plane. Each pole will be associated with one pseudomode. Reverting to Schrodinger picture amplitudes 
via ci{t) — ci(t)e^"^i* etc., pseudomode amplitudes associated with each pole of D are introduced as defined by: 

t 

izit' 



bi{t) = -il7V=Mie-'^'* / dt'e'^'' ci{t') (10) 



Prom the definition of hi and by substituting the form (Eq.(^)) for the kernel G that involves the poles oi D, it is 
not difficult to show that the excited atomic amplitude and the pseudo mode amplitudes satisfy the following coupled 
equations: 



dt 

, dbi{t) 
* dt 



^ zibi{t) + Kiciit) (11) 



where /C; = Cl^—in are pseudomode coupling constants. In general, the residues r; are not pure imaginary, so the 
pseudomode coupling constants are not real. 

The important point is that the atom plus pseudomodes system now satisfies Markovian equations (Eq.(pT|)). With 
a finite (or countable) set of pseudomodes, the original atom plus structured continuum has now been replaced by a 
simpler system which still enables an exact description of the atomic behaviour to be obtained. Exact master equations 



involving the pseudomodes have been derived, and in general the pseudomodes can be coupled (see 51 1). There are 
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however, difficulties in cases where the pseudomode couphng constants are not real, which can occur in certain cases 
where there are several pseudomodes. Apart from the general difficulty associated with situations where (apart from 
having simple poles) the reservoir structure function D is not analytic in the lower half complex plane, there is 
also considerable difficulty in extending the above theory to treat cases where multiple excitations of the structured 
reservoir occur, such as when the two-level atom is replaced by a three-level atom in a cascade configuration with the 
atom is initially in the topmost state. The problem is that applying the usual essential states approach leads to two 
(or moreVphoton states now appearing in the state vector (cf. Eq.(0)), and the resulting coupled amplitude equations 
(cf. Eq.(y)) do not apear to facilitate the sucessive formal elimination of the one, two, .. photon amplitudes, as is 
possible in the single photon excitation case treated above. It is therefore not clear how pseudomode amplitudes can 
be introduced, along the lines of Eq.(|l0|) so the pseudomode method has not yet been generalised from its original 
formulation to allow for multiple reservoir excitations. 

III. FANO DIAGONALISATION FOR A QUASI MODE SYSTEM 

A. Description of the approach 

The case of multiple excitation of a structured reservoir involves systems more complex than the two level atom 
treated above. It will be sufficient for the purpose of linking the pseudomode and Fano diagonalisation methods to 
consider single multilevel atomic systems, although multiatom systems would also be suitable as both systems could 
result in multiphoton excitations of the quantum EM field. Accordingly the two level Hamiltonian given as the second 
term in Eq.(|l|) is now replaced by the multilevel atomic Hamiltonian: 

Ha = Y^ Vkhujk {(rt^^k - ^k(^k) (12) 
fc 

The index k represents an atomic transition associated with a pair of energy levels (fc = {u, I}) with energy difference 
hujk- The quantities rjk are numbers chosen so that Ha equals the atomic Hamiltonian, apart from an additive 
constant energy. For example, in a two level atom f] = \ for the single transition, whilst in a three level atom in a V 
configuration with degenerate upper levels ?7i = ?72 = for the two optical frequency transitions, and ??3 = for the 
zero frequency transition. Details are set out in Appendix The atomic transition operators are = = (o"^)^- 

As the Hamiltonians for other fermionic systems can also be written in the same form as in Eq. (|2|) , the treatment is 
not just restricted to single multilevel atom systems. 

As indicated in Section an important pseudo mode situation is where the reservoir structure is due to the 
presence of a discrete, system of (quasi) modes which are coupled to other continuum (quasi) modes with slowly 
varying coupling constants. This important case applies to atomic systems coupled to the quantum EM field in high 
Q resonant cavities, such as microspheres or microcavities. The Fano diagonalisation method is then based around 
the idea that the structured reservoir of quantum EM field modes can be described in two different ways, which will 
now be outlined. Figure ^ illustrates these two descriptions, along with that involving pseudo modes. 

1. Quasi modes 

The first approach is to treat the quantum EM field in terms of a quasi mode description ||5^,|5^ . The quasi modes 
behave as coupled quantum harmonic oscillators. These are to consist of two types, the first being a set of discrete 
quasi modes, the second being a set of continuum quasi modes. In a typical structured reservoir situation for the area 
of cavity QED the quasi modes represent a realistic description of the physical system. The discrete modes would 
be cavity quasi modes — one for each cavity resonance and appropriate for describing the EM field inside the cavity, 
the continuum modes would be external quasi modes which are describe the field outside the cavity. The interaction 
between the discrete and continuum quasi modes will be treated in the rotating-wave approximation assuming slowly 

6l| ]. Rotating-wave approximation couplings between the discrete quasi modes are 
also included, but couplings between the continuum quasi modes are not included — such couplings can be removed 
by pre-diagonalisation. For the quasi mode description the field Hamiltonian is given by: 

Hp = fiujia\ai + hvijolaj 



varying coupling constants [pO 57 
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nWii A)alb{A) +h.c. 

+ y dApe(A)nASt(A)S(A) (13) 

where Oi, a| are the annihilation, creation operators for the discrete quasi mode i, uJi is its frequency, ^(A), {A) are 
the annihilation, creation operators for the continuum quasi mode of frequency A, the coupling constants between the 
i,j discrete quasi modes are Vij {vij = Vji), whilst the quantity Wi(A) is the coupling constant between the i discrete 
and A continuum quasi modes. The integrals over the quasi continuum frequency A involve a quasi continuum mode 
density pc{A). Both l/Fi(A) and Pc(A) are usually slowly varying. The discrete quasi mode annihilation, creation 
operators satisfy Kronecker delta commutation rules, whilst those for the continuum quasi mode operators satisfy 
Dirac delta function commutation rules: 

[flj, a]] = Sij 

[S(A),6t(A')]-^(A-A')/pc(A). (14) 

The pc factor on the right hand side gives annihilation and creation operators which are dimensionless. 

For the quasi mode description the interaction between the atomic system and the quantum EM field will be given 
in the rotating-wave approximation and only involve coupling to the discrete quasi modes. This would apply for the 
typical structured reservoir situation for the area of cavity QED in the familiar case where the atoms are located 
inside the cavity. The energy of an excited atom escapes to the external region in a two step process: first, a photon 
is created in a discrete (cavity) quasi mode via the atom-discrete quasi mode interaction, then second, this photon 
is destroyed and a photon is created in a continuum (external) quasi mode via the discrete-continuum quasi mode 
coupling. For the quasi mode description the atom-field interaction will be given as: 

^AF^Y^Y. i^^i^'^^^i + ii-c-) (15) 

k i 

where Xki is the couphng constant for the k atomic transition and the i quasi mode. 



+ ^ / dAp,(A) 



2. True modes 

The second way of describing the quantum EM field is in terms of its true modes [58|j59|] . The true modes behave 
as uncoupled quantum harmonic oscillators. These modes are also used in cavity QED and are often referred to as 
"universe modes" . The pseudomode theory presented in Section || is also based on true modes. For frequencies near 
to the cavity resonances these modes are large inside the cavity and small outside, for frequencies away from resonance 
the opposite applies. The distinction between true modes and quasi modes is discussed in some detail in recent papers 
|p6,^ and their detailed forms and features in the specific case of a planar Fabry-Perot cavity are demonstrated in 
pO | In terms of true modes the field Hamiltonian is now given in the alternative form as: 

Hp = j dio p{u))nijjA\ijj)A{ijj) (16) 

where A{uj), {'jj) are the annihilation, creation operators for the continuum true mode of frequency w. The integrals 
over the quasi continuum frequency lo involve the true continuum mode density p{oj), which is not in general the same 
function as pc{A). It is also not necessarily a slowly varying function of uj. The continumn true mode annihilation, 
creation operators satisfy Dirac delta function commutation rules: 

[A{oj),A\uj')]=6{uj-J)/p{uj) (17) 

In all these Hamiltonians the coupling constants have dimensions of frequency, whilst the annihilation and creation 
operators are dimensionless, as are the atomic transition operators. 



3. Relating quasi and true modes 



As will be demonstrated in Section IIIB, Fano diagonalisation involves determining the relationship between the 
true mode annihilation operators A(uj) and the quasi mode annihilation operators di and b{A). The A{uj) will be 
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written as a linear combination of the ai (sum over i) and 6(A) (integral over A) (see Eq.(p^) below), which involves 
the functions ai(ijj) and P{uj,A). This relationship can be inverted to give the as an integral over io of the A{uj) 
(see Eq.(^2|) below). This enables the true mode form of the atom- field interaction to be given as: 

Haf = J2J2 f"^ P^"") (^^^^"* + h.c.) (18) 

k i '' 

Comparing Eqs.(p^ and ( p^ we see that the atom-true mode coupling constant g''{uj) (for the k atomic transition 
and the lu true mode) is given by the expression: 

/(c^) =^Afc,a,(cc>). (19) 

i 

This can be a complicated function of lu in a structured reservoir, as will be seen from the forms obtained for the 
function ai{uj) (for example, Eq.(|67|)). This expression for the atom-true mode coupling constant is one of the key 
results in our theory, and enables the pseudomode and quasi mode descriptions of decay processes for structured 
reservoirs to be related. Note that the true mode coupling constant now involves two factors: the atom-quasi mode 
coupling constant Xki and the function ai{uj) that arises from the Fano diagonalisation process. 

For the situation where only a single atomic transition k is involved, the equivalent reservoir structure function 
would be given by: 

D'^{u;)^Cp{iu)\g\cj)\' (20) 

where C is the normalising constant, which for convenience we will set equal to unity as it does not contain any uj 
dependence. This expression will be used to compare the results from the quasi mode approach to those of the present 
single quantum excitation pseudomode theory. As we will see, the true mode density cancels out. 

Finally, athough our results are still correct for cases where the quasi mode density /9c(A) and the coupling constants 
Wi(A) are not restricted to being slowly varying functions of A, their utility where this is not the case is somewhat 
limited. The theory is mainly intended to apply to the important pseudo mode situation where the reservoir structure 
is actually due to the presence of a discrete, system of quasi modes which are coupled to other continuum quasi 
modes via slowly varying coupling constants. For example, the quantum EM field in high Q resonant cavities can 
be accurately described in terms of the quasi mode model which has these features, the discrete quasi modes being 
the cavity quasi modes (linked to the cavity resonances) with which the atoms inside the cavity interact, and the 
continuum quasi modes being the external modes. 

As pointed out previously, the structured reservoir can be any set of bosonic oscillators, not just the quantum EM 
field. The above treatment would thus apply more generally, and we would then refer to discrete quasi oscillators, 
continuum quasi oscillators or true oscillators. The physical basis for a quasi mode description of the reservoir of 
bosonic oscillators will depend on the particular situation; in general they will be idealised approximate versions of 
the true modes. 



B. Diagonalisation of the quasi mode Hamiltonian: dressing the quasi mode operators 

1. Basic equations for Fano diagonalisation 

We start with a multiple quasi-mode description of the quantum EM field, for which the Hamiltonian is given above 
as Eq. ([l3|) . This Hamiltonian can also be written in terms of the true mode description as in Eq. (|l6|), and the 
problem is to relate the true mode annihilation operators A{lo) in terms of the quasi mode annihilation operators 
and 6(A). In view of the rotatin g w ave approximation form of the Hamiltonian, the quasi mode creation operators are 
not involved in the relationship |5q ] . Fano diagonalisation for the non-rotating wave approximation has been treated 
for the case of a single mode coupled to a reservoir in Ref. In making a Fano diagonalization we will follow 

the lines of Ref. jlj] (Section 6.6 on dressed operators), rather than Ref. [Q, but note that a new feature here is the 
presence of the mode-mode coupling term in the Hamiltonian Eq. (p^) . In addition, we explicitly include the mode 
densities from the beginning. The physical realisation of the quasi mode model for the EM field really determines 
the quasi continuum mode density Pc(A), just as it does the coupling constants Vij,Wi{A) and Xki- It is therefore 
important to be able to find the Pc{^) dependence of quantities such as the reservoir structure function D{u>) (as we 
will see, the final expression (Eq.(pO[)) for the latter does not involve the true mode density /^(a;)). It is of course 
possible to scale all the other quantities to make p = pc = 1, and then rescale afterwards to allow for the actual 
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p, Pc that apply for the system of interest, but this would lead to much duplication of the results we present. For 
completeness, the scaling is set out in Appendix 

From the form of the true mode Hamiltonian in Eq. (^6|) and the commutation rules Eq. (|l^) to be satisfied by the 
A{lu), it is clear that the true mode annihilation operators are eigenoperators of the quantum field Hamiltonian Hp 
and must satisfy: 



(21) 



In general, the true mode annihilation operators A{ll!) can be expressed as linear combinations of the quasi mode 
annihilation operators a, and 6(A) in the form ( [p7|,p6[): 



AiLo) 



(22) 



where ai{uj) and P{lu,A) are functions to be determined, and which are dimensionless. This form for A{lu) is 
then substituted into Eq.(21) and the commutator evaluated using the quasi mode form Eq.(p^) for Hp and the 
commutation rules in Eq.(L4). The coefficients of the the operators at and 5(A) on both sides of Eq. (pT|) are then 
equated, giving a set of coupled equations for the ai{uj) and l3{u!, A). These are: 



[ui - uj)a^{uj) + + / dApc(A)/?(tJ, A)W*{A) = 

(A - u)I3{lo, A) + ^ W,{^)a,{Lo) = 



(23) 
(24) 



To solve Eqs.(p3[), ( |24| ) for the unknown ai{u}) and (3{ijJ, A) wc first solve for (3 in terms of the a^. This gives: 

V 



Pico, A) 



-A 



+ z{lu)S{uj - A) 



J2w,{A)a,iLo), 



(25) 



where z{uj) is a dimensionless function yet to be determined. This expression is then substituted into Eq.(|23|) to 
obtain a set of linear homogeneous equations for the at (w) in the form: 



{uji - uj)ai{uj) + ^ Vjiaj{uj) + ^i^yaj(w) + ^ W*{uj)W.j{LL>)pc{uj)z{uj)aj{uj) = 0. 

j^^i 3 3 



(26) 



In these equations, a frequency shift matrix Fij{uj) appears, which involves a principal integral of products of the 
discrete-continuum quasi mode coupling constants together with the quasi continuum mode density. This is defined 
by: 



W*{A)W,{A) 
uj~ A 



and satisfies the Hermiticity condition Fji = F*^ . 
Equation ( p6|) can be written in the matrix form 

mo; = 

where the column matrix a. = {ai(a;), q;2(w), a^{uj), ...}^ and the square matrix m is given by: 
my (w) = {uji ~ Lo)6ij + (1 - 5^j)vJl + Fij{uj) + W* {uj)Wj{Lo)pc{io)z{uj). 



(27) 

(28) 
(29) 



2. Solution of equations for amplitudes Qi(aj) and /3(a;, A) 



The approach used to solve these equations is as follows. It is clear that Eq. (p8[) can give an (unnormalized) solution 
for OL in terms of the function pc{ijj)z{uj). We can now use Eq. (p8|) itself to obtain the expression for pc{u>)z(u)), subject 
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to the assumption that the quantity Wi(w)ai(w) is non-zero. This assumption will be verified a posteriori from 
the normalisation condition for the ai{uj), which will follow (see below) from the requirement that the form for the 
A{ijj) given in Eq.(|2^), satisfies the commutator relation_L4(w), — 5{uj — uj') / p{uj), [Eq.(|l7|)]. This indeed leads 

to a non-zero expression for Wi{uj)ai{uj) ^ (see Eq.(p9[) below). After finding both pc{u})z{u}) and Wi{bj)ai{u}) 
the results can be substituted back into the equations (p6D. By eliminating the factor Wi{ijj)ai{u)) from the last 



term in Eqs.(26), we obtain a set of inhomogeneous linear equations for the ai{uo), which can then be solved for the 
ai{Lo) (and hence /3(w, A)). 

The general expression for Pc{^) z{Ld) can be obtained from the matrix equation (pq). With E the unit matrix we 



introduce the square matrix Jl, the column matrix W* and the row matrix W"^ via: 

r2jj(cj) = uJiSij + Sij)vy^ + Fij{uj) (30) 



and W*(w) = {W^{uj),W:^{uj),W^{uj),...Y , W^(w) = {Wi{uj),W2{uj),W-i{u)), ...} , and then write Eq.(|8l) in the 
form: 

{-{ujB-n) + pc{uj)z{uj)W*W'^)oL = Q. (31) 

Now the matrix Jl is Hermitean and positive definite, having real eigenvalues close to the real and positive uji. The 
matrix — Vl can be hence assumed to be invertible, so by multiplying Eq.(^l|) from the left by W-'"(wE — Jl)^^ we 
see that: 

(-1 + Pc(w)z(w)J(w))W^a = (32) 

where the function J(ui) is defined by 

7(0-) = W^(wE- r2)"^W*. (33) 

Now the quantity 'W^ a. is equal to Wi{ijj)ai{Lo), which is assumed to be non zero for reasons explained above. 
This means that (—1 -I- Pc{lo)z{lj)J{ijj)) = 0, and this gives for Pc{lo)z{uj) the general result: 

Pc{io)z{u)=\^WM{^^-^{^)%^W;{u)^ , (34) 

which only involves the various coupling constants and angular frequencies, along with the quasi continuum mode 
density. In general the w dependence of the result for pc{uj)z{uj) is complicated, since both the coupling constants Wi 
and the matrix (via the matrix F) will depend on lj. In some important cases however, their uj dependence can be 
ignored. 

As indicated previously, Eqs.(|2^) or ( |28| ) only determine the ai{uo) (and hence P{lo, A)) to within an arbitary scaling 
factor, as can be seen from their linear form. The normalisation of the solutions is fixed by noting that we need A{uj), 
Eq. (p2[) , to satisfy the commutator relation [j4(ti;), A'''(w')] = 5(bj — u}')/ p{lo), [Eq.(p^)]. This leads to the condition: 



^ a,{uo)a*{J) + j dAp,iA)f3iu;, A)/3*(c^', A) = 5(c^ - ij')/piu;). (35) 
Then substituting for (3{lu,A) from Eq.(p5|) and using Eq. (p7|) , we find after considerable algebra that 

V 



UJ — UJ' 

I] 



+z*{uj')p,{uj')W,{uj')W;{uj') - z{uj)p,{oj)W,{uj)W;{oj)\ = 5{uj - oj')/p{uj). (36) 

Note that we have used certain properties of the principal parts and delta functions (sec, for example Ref. []l2| ): 

5{uj - A)5{uj' - A) = 5{uj - uj')5{uj - A) = 5{uj - uj')5{uj' - A) 

V V 

-6(uj-A)^—^ b{uj-A) 



w' — A uj' ~ UJ 

V V _ V ( V v_ 

uj-A uj' - A ^ UJ - uj' \uj' - A uj-A 



+ n^S{u! - A)S{uj' - A) (37) 
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to obtain the last equation. We then also use 

r 



- ■ {uj-uj')^l (38) 



UJ ~ UJ' 

along with Eq.(p6|) to substitute for F*^{ijj')a*{Lo') and Fji{ijj)ai{ijj) and obtain finally: 

2 

(39) 



(w)ai(w)a* (w) 



Y,W^{uJ)a^{uJ) 



p(c.)p,(u;)(7r2 + |z(o.)|2)- 



This fixes, albeit with the coefficients Wi{Lu)^ the normalization of the ai(aj). Note the appearance of both mode 
densities in the result. Finally, with a suitable choice of the overall phase we can fix the result for the important 
quantity Wi[uj)ai[uj) to be: 

E W.{u,)<^.{u^) = / t . t\i ^ - I (40) 



Having obtained this result for '^^Wi{uj)ai{Lij) we then substitute back into the equations (^6|), eliminating this 
factor from the last term to give a set of inhomogeneous linear equations for the Q!i(aj): 

[uj ~ UJ^)a^[uJ) -^Vj,aj{uj) - 2_^F,jaj{uj) = ■ — . (41) 

j VP(^^)Pc(t^)(7r + «z(w)) 

After some algebra, introducing the matrix fl^uj) from Eq.(|30|) and then substituting from Eq.(|3^) for {pc{uj)z{ll!))~^ , 
the last equations can be solved for the ai(Lu), giving the solution in matrix form as: 



aiuj)= , ,„,^, — -— — — --(a)E- n(tj))-iW*(tj). (42) 

^ ' Y p{uj) (l-i7rpc(t^)W^(w)(wE-f2(w))-iW*(w))^ ^ " ^ ' ^ ' 

In this result all the terms that in general depend on lu are explicitly identified. It is also convenient to write the 
inverse matrix in terms of its determinant and the adjugate matrix via: 

(wE - n{Lu))-^ = (wE - r2(w))^'°'V|t^E - n{uj)\ (43) 

and then the solution for a{uj) becomes: 



f N . /PC(W) 1 

a.[uj) = —I ' 



p{uj) (|t^E - n{uj)\ - i7rpe(w)W'r(w)(wE - f1(w))^^^W*(u;)) 
x(cjE-0(cj))^^''W*(w). (44) 

The result for the expansion coefficient f3{uj,A) then follows from Eq.(p5|) and substituting for Pc{uj)z(lu) from 
Eg . (^4|) . After some algebra we find that: 

1 [Siiu^ A) + ^p^{uj)W^{u;){ij-E~n{u;))-^W*{ij)] 
" ~Vp(^)Pc(^) (1 - i7Tp,{Lj)WT{u;){cu-E - n{u;))-^W*{u;)) ^^^^ 

We see that the solutions for the ai{uj) and /3(w, A) only involve the various coupling constants and the mode densities. 

3. Coupling constants and reservoir structure function 

Introducing the column matrix Afc = {A^i, Afe2, A^a, ...}^ the expression ( |T9| ) for the coupling constant g^{ijj) can 
be written as: 
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1 



Y p{u;) (|tjE-n(w)| -i7rpc(w)W^(w)(wE-r2(w))'4^-^W*(w)) 
X A^wE - 0(cj))^^''W*(cj) 

*Y p(c^) P„(w) 
where the functions Pn{io) and Qjj_]^(ti>) are defined by: 



(46) 
(47) 



ij 



(48) 



(49) 



In the case where the uj dependence of the quantities pdco) , Fij (uj) and Wi{u!) can be ignored, P„(aj) and Qn-ii^^) 
would be polynomials in co of degrees n and n — 1 respectively, as will be seen in Section [V . 
The reservoir structure function can then be expressed as (C = 1): 



(50) 



where we note the cancellation of the true mode density p(a;) and the proportionality to the quasi continuum mode 
density Pc(w)- The significance of the p{uj) cancellation will be discussed in Section IIIC. There is however further 
dependence on the quasi continuum mode density within the function P„(w), as can be seen from Eq.ji^). The role of 
this dependence will be discussed in Section IV when we have obtained expressions for the reservoir structure function 
for specific cases. 

To sum up: if we are given the Hamiltonian in the quasi mode form Eq.(O), we can obtain the true mode operators 
( p^ ) which satisfy the eigenoperator condition Eq.(|2^). The coefficients ajuj) are found by solving ma — 0, Eq. (p^); 
the function z(aj) occuring in m is obtained from Eq.(|28|) and given by Eq.(^^. The solutions for ajjui) are scaled 
in accordance with Eq.(|35|) and the normalisation for the quantity Wi(w)ai(w) is given in Eqs.(|39|), (pO|). The 
normalised solutions for Q!i(Lj) are obtained as Eqs.(42) or (^). The coefficients P{u), A) are then found from Eq.([25|) 
and the result given in Eq.(Ea). The true mode coupling constant g^{Lo) and the reservoir structure function D^{Ld) 
are obtained as Eqs.(^) and (|50|). These results involve the functions Pnioj) and Q'^_i{lo) defined in Eqs.(48) and 
(p9|). The results depend on the quasi continuum mode density pc as well as on the various coupling constants and 
angular frequencies. It should be noted that a unique expression has been obtained for z{ijj), and hence for the aiiuj) 
and (3{uj, A), even though the determinental equation |m| =0 might appear to give anything up to n solutions, where 
n is the number of discrete quasi modes. This feature is due to the specific form of the matrix m that is involved. 
The overall process amounts to a diagonalization because the EM field Hamiltonian in the non-diagonal quasi mode 
form is now replaced by the diagonal true mode form given by Eq.(16). 



C. Inverse diagonalization: undressing the true mode operators 

We can also proceed in the opposite direction from Fano diagonalization: that is, we can also find the quasi mode 
operators and 6(A) in terms of the true mode operators A{uj). In general ( [ |57| , p6t ) the quasi mode annihilation 
operators and b{A) can also be expressed as linear combinations of the true mode annihilation operators A(uj) in 
the form: 

di = j duj p{uj)^i{uj)A{uj) 

b{A)^ J dLup{Lu)S{A,uj)A{uj) (51) 

where the functions 7i(i^) and 5{A,lu) have to be determined. These can be obtained in terms of the ai{uj) and 
/3{lu, A) by evaluating the commutators [A(a;), d|] and [^(w), 6(A)t] using the basic commutation rules Eqs.(p7|), (pH). 
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For the first commutator: on substituting for A{u;) from Eq.(22) we obtain ai{uj), on the other hand, substituting 
instead for di from Eq.(|5l|) gives j*{uj), and hence = 7*. Carrying out a similar process for the second commutator 
gives the result [3 = 5* and thus: 



duj p{Lj)a* {lo)A{ljj) 

b{A) = I dLjp{uj)(3*{A,uj)A{uj) (52) 



As has been already described in Section pi A| , the first of these two equations enables us to relate the two descriptions 
of the atom-field interaction given in Eqs.(|15|) and (^8|). Ultimately, the key expression we have obtained in Eq.(p^) 
for the atom-true mode coupling constant rests on this result. As we will see in Section this enables us to relate 
pscudomodes to the discrete quasi modes. 

As a final check of the detailed expressions, in Appendix ^ we start with the field Hamiltonian in the quasi mode 
form Eq. (p^) , then substitute our solutions for ai{uj) and f3{Lj, A) into the expressions for d^ and b(A) given in Eqs.(^2[). 
On evaluating the result, the Hamiltoni an in the true mode form Eq.(^6|) is obtained — as required for consistency. 

It has already been noted in Section IIIE| that the final expression for the reservoir structure function D^{uj) in 



terms of quasimode quantities is independent of the true mode density p(w). Also, we have found no equation that 
actually gives an expression for p{lo) in terms of the quasi mode quantities, including the continuum quasi mode 
density Pc{A) - a somewhat surprising result. The true mode density therefore does not play an important role in 
the quasimode theory. The reason for this is not that hard to find, however. The theory can be recast with both 
the p{uj) and Pc{A) factors incorporated into the various operators and couphng constants. In Appendix ^we show 
that p{u>) and Pc{A) can be scaled away to unity. For example, from Eqs.(^2|), ( ^ ) and (|2|) we see that the true 
mode annihilation operator is proportional to 1/ yj p(w), the other (operator) factor only depending on quasi mode 
quantities. Hence (as in Appendix ^ we may scale away the p(uj) dependence via the substitution: 

where A*^*) (w) is independent of p. If this substitution is made then the field Hamiltonian is given by: 

Hf = J dcj (^)^('i) (^) (54) 

without any p{lo) term. 

IV. APPLICATIONS 

A. Case of a single quasi mode 

For this case no coupling constant between discrete quasi modes is present and we may easily allow for a non zero shift 
matrix element Fu and for non constant WjJA). Noting that {luE — Ct{Lj))^^ = 1 and jojE — r2(aj)| = uj — uji — Fii{uj), 
a simple evaluation of Eqs.(|3^), ( ^ ) and (|4^ ) gives the following results: 

Pc(c.).(c.) = (55) 



g'^ito) = Afciai(w) 



y p{lo) lo - LOl- Fii{uj) - iiT pc{oj)\Wi{ijj)\'^ 
In terms of a frequency shift Awi and half-width ^ defined as: 
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Awi(u;) = Fii{lu) 



(57) 



^ = np,iu;)\W,iLo)\^ (58) 

the reservoir structure function (see Eq.(^o|)) for the situation where only a single atomic transition k is involved is 
then found to be (C = 1): 

pkr. \Xki\'-r{u;)/2n 

In the situation where the quasi mode density Pc(A) and the coupling constant Wi(A) are slowly varying functions 
of A, these quantities can be approximated as constants in the expressions for the frequency shift and width. The 
reservoir structure function is then a Lorentzian shape with a single pole in the lower-half plane at loi + Awi — iT /2 
corresponding to a single pseudomode. Thus the single discrete quasi mode is associated with a single pseudomode, 
whose position zi is given by u)i + AtJi — iV /2 in terms of quasi mode quantities. 

B. Case of zero discrete quasi mode-quasi mode coupling and flat reservoir coupling constants 

The theory becomes rather simpler if there is no coupling between the discrete quasi modes, that is: 

=^ (60) 

This could be in fact arranged by pre-diagonalising the part of the Hamiltonian Hp that only involves the discrete 
quasi mode operators. Thus we write 

huJialcLi + ^ hvija\dj (61) 

in the form 

i 

via the transformation 

c,^^U^jaj. (63) 
j 

where U is unitary. The last equation can be inverted to give the di in terms of the Ci and the result substituted in 
other parts of Hp (Eq.(p^) and Haf (Eq.(p^). The original coupling constants Xki and Wi{A) would be replaced 
by new coupling constants via suitable linear combinations involving the matrix U, and these generally would have 
similar properties (e.g., flatness) as the original ones. 

The idea of replacing the structured reservoir of true modes by quasi modes, in which the continuum quasi modes 
constitute a flat reservoir, implies that the discrete-continuum quasi mode coupling constants Wi{A) and the quasi 
continuum mode density /3c(A) are slowly varying functions of A. This results in the shift matrix Fij elements being 
small, so it would be appropriate to examine the case where they are ignored, that is: 

F,, =^ (64) 

with both pc and the Wi are assumed constant. 

For the case vtj = 0,Fy = 0,pc{A) = pc and Wi{A) = Wi (constants) the quantities involved in the inverse of the 
matrix cuE — n{Lu) are: 

jwE — Jl(a;)| = (ll) — uji){lo — lu2)---{'-^ ~ ^n) 
(wE - r2(w))^-°"' = (tj - uji){uj - tj2).-.(w - uji-i){uj - Wi+i)...(tj - w„)(5y . (65) 



A straightforward application of Eqs.(34) and (M) leads to the simple results 
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^=-H = |Ei^| (66) 

^^(^) = -^,^^W: " ~ ~ "r:^^" ' ~ (67) 

where the function P„([j), (which is defined in Eq.(^8|)) is now a polynomial of degree n, whose roots are designated 
as ^i. It is now given by: 

P„(w) = (UJ - LlIi){uJ - W2)---(w - iOn) 

-inpc'^ \Wj\^{uJ ~ LOl)...[uj - tJj_l)(w - UOj + l)...(uj - LOn) 
3 

= (cc>-a)(^-6)-(^-en). (68) 

For the true mode coupling constants g^{Lo), the general result in Eq.(^7|) can be applied to give: 



p{uj) {Ul - Ci)(w - £,2)--{l^ - C«) 



(69) 



where the function Q^^_i(uj) (which is defined in Eq.(^)) is now a polynomial of order n—1, whose roots are designated 
as Oi. It is now given by: 

Qn_l{uj) = ^ XkiW*{uj - UJl){uJ ~ UJ2)...{UJ - UJi-l){u! - UJi+i)...{uJ - UJn) (70) 

i 

= Skiu; - 9i){uj - 62). ..{u - 0„_i), (71) 
where 5^ is a strength factor defined as: 

^fc=^Afc,W,* (72) 



The reservoir structure function D [u) (see Eq.(50)) for the k transition is then given by (C = 1): 



r,k., ^ ic, |2 l(^-gl)(^-^2)...(c^-gn-l)P 

|(tJ-$l)(cj-^2).-.(t^-Cn)r 

Since products of the form {uj — ^ C*) can be written as (w — Ref)^ + (Imi^)^, the behaviour of the reservoir 
structure function D'^ = p{uj)\g''{uj)\'^ (see Eq.(pO|)) as a function of to is now seen to be determined by the product 
of n Lorentzian functions associated with |P„(a;)P with the modulus squared of the polynomial of degree n—1 given 
by |Q^_2(<x')p. The quasi continuum mode density merely provides an uninteresting multiplicative constant, except 
insofar as it is involved in expressions for the width and shift factors. In the case where there are n discrete quasi 
modes, then irrespective of the location of the roots of the polynomial equation P„(ijj) = 0, the reservoir structure 
function D^{lu) for a single quantum excitation has n poles in the lower half plane, each corresponding to either or 
As there are n roots when n discrete quasimodes are present, we see that each discrete quasi mode corresponds 
to one of the n pseudomodes, whose position Zi is equal to or to Thus, for the case here where the coupling 
constants and the quasi continuum mode density are independent of frequency, the feature that leads to a pseudomode 
is the presence of a discrete quasi mode. 

C. Case of two discrete quasi modes 

The results in the previous subsection can be conveniently illustrated for the case of two discrete quasi modes. For 
simplicity we will again restrict the treatment to the situation where vi2 = 0, i^ij = 0,pc{A) = pc and VFi(A) = Wi 
(constants), and just consider a two-level atom, so only two coupling constants Ai, A2 are involved. In this case the 



atom-true mode coupling constant can be obtained from Eq.(69) and is: 
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g(uj) = -i.—— — — (74) 

V p(w) (u;-^i)(tj-C2) 

where and the roots ^1.2 of P2{lo) — are given by: 

{XiW* + A2 ) {XiW^ + A2 Wi ' 

and 

ei,2 = \{{uoi + + i-^PciWxf + 11^2!')} 

±iV{K - ^2) + - W2Y)Y - ^^^pI\W^^\W2Y (76) 

It will also be useful to introduce widths Vi defined by: 

r, = 27rp,|V7,p (77) 

and which can be later identified (see Section |^ as the discrete quasi mode decay rates (Eq.(|9S|)). These results will 
be now examined for special subcases. 

1. Special subcase: Equal quasi mode frequencies 

In this case we choose: 

wi = a;2 = (78) 

and find that: 

U)Q = LOc 

6.2 =c^c,wc + «7rp,(|VFi|2 + I VF2n (79) 

giving for the atom-true mode coupling constant: 

— {XiWl + A2Ty2*) 
p{uj) L0-LOC + t^pcdmp + IVF2P) '■ 

and for the reservoir structure function: 



g(^) = -h/ ,. , /iw |2^lw/J2^ ' (80) 



U{LJ) - Pc- ^2 I a-p 1 1 /o^2 ^^'^ 



This corresponds to a single pole in the lower half plane for the reservoir structure function (see Eq.(|2C|)) and thus 
only results in a single pseudomode, albeit for a case of two degenerate discrete quasi modes. 

2. Special subcase: Equal quasi mode reservoir coupling constants 

In this case we choose: 

Wi^W2 = W (82) 



and find that: 



2''^! + ^2) 
(Ai - A2) , , 
^"^^ 2(Ai + A2) ^"^'"^^ 

6.2 = + U2) + t7Tp,\W\^ ± ^^{UJI - LJ2Y - 4^2p2|Vl/|4 (83) 
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Here ujq has been written in terms of the quasi modes centre frequency ujc and a frequency shift Au>c depending 
on the difference between the two atom-discrete quasi modes couphng constants Xi and the discrete quasi modes 
detuning. There are now two regimes depending on the relative size of the discrete quasi modes separation jwi — W2I 
compared to the square root of the quasi continuum mode density .^/p^ times the reservoir couphng constant W. 
Equivalently, the regimes depend on the relative size of the separation jcji — lu2\ compared to the width factor (decay 
rate) T = Ti = = 27rpc\W\^. 

a. Regime 1: Large separation \uji — uj2\ > F Adopting the convention that cui < UJ2 we can write 

^^{uji - LJ2Y - Att^pI\W\^ = ^(^2 - c^i) - Alor (84) 
where I^lur is a reservoir induced frequency shift. The atom-true mode coupling constant now becomes: 

( ^_■ (Ai + \2)W*{UJ -UJC- A^c) , . 

and the reservoir structure function is then: 

^_ |Ai + A2p(r/27r)(c^-^c-A^c)2 

^ ^ ((c^-cj2 + Ac^fl)2+r2/4)((c^-c^i-Ac^;^)2 + r2/4) ^ ' 

The reservoir structure function D (see Eq.(p6|)) will be zero at the shifted centre frequency uJc + Ai^c- There are 
two poles in the lower half plane leading to Lorentzian factors centred at frequencies UJ2 — Atj/j and uJi + Aujr and 
which have equal widths 2npc\W\'^ . We note that the effect of the coupling to the reservoir is to decrease the effective 
discrete quasi modes separation by 2Aujr. 

h. Regime 2: Small separation \loi — uj2\ < T We now write 

i V4^2^2|p^|4 _ _ ^2)2 - ^Pc\W\''{\ - A/r) (87) 

where A/p is a fractional change in width factors associated with discrete quasi mode separation. The atom-true 
mode coupling constant now becomes: 

. ^ rpT (Ai + A2)W^*(a>-^c- Ac^c) 

^^''> '\j p{Lo){uo-ujc-'2iTrpc\W\m-\Afr)){uo-uc-iTipc\W\^Afry ^ ' 

and the reservoir structure function is: 

j^, , \\^ + X2?{T/2i,){uj-ujc-AoJc? , . 

The reservoir structure function D (see Eq. (p9[) ) will again be zero at the shifted centre frequency wc' + Acjc- There 
are two poles in the lower half plane leading to Lorentzian factors both centred at the same frequency wc, but which 
have unequal widths 27rpc|M^P(l — ■^A/r) and 7rpc|M^pA/r. If A/r ^ 1 one width is much smaller than the other. 

In their work on superradiance in a photonic band gap material Bay et al |]65| assume as a model for the mode 
density a so-called Fano profile of the form: 

with the two-level atom coupling constant g{uj) given by a slowly varying function proportional to It is interesting 
to note that the reservoir structure function related to their theory is of the same form as that obtained here from 
Eq.(^9|) if the following identifications arc made: 

q — > -Aloc 
\k ^ 2iTp,\WWl - \Afr) 

\p ^ vrp.lM^pA/r (91) 

For situations such as atomic systems coupled to the field in high Q cavities the physics is different of course, with 
the resonant behaviour in the reservoir structure function being due to the atom-true mode coupling constants rather 
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than the reservoir mode density (which we assume is slowly varying). Nevertheless, our two discrete quasi mode 
model — with equal reservoir coupling constants W that are large compared to the discrete quasi modes detuning 
1^1 — uj2\ — does provide an equivalent physical model for the photonic band gap case that Bay et al treated, the lack 
of which was commented on in the review by Lambropoulos et al ||2l|] . 

The band gap case was also treated as a specific example by Garraway in the original pscudomodc theory 
paper. A model for the reservoir structure function was assumed in the form of a dilference between two Lorentzians: 

Fi r2 

where the weights wi , W2 satisfy wi — W2 — I. Again, apart from an overall proportionality constant this same form 
can be obtained here (see Eq.(^)) for the reservoir structure function D if we choose the atom-discrete quasi mode 
coupling constants Ai, A2 to be equal (so that the frequency shift Auc is zero): 

Al — A2 

Aujc = (93) 

and where the following identifications are made: 

^T2^7rp,\WfAfr 

Wi 



1 - 



ir 



1-A/i 
^A/r 



r 



- (94) 

As will be seen in the Section ^ the existence of unusual forms of the reservoir structure function (such as the 
presence of Lorentzians with negative weights) does not rule out Markovian master equations applying to the atom- 
discrete quasi modes system. Thus, for the situation of single quantum excitation, where the pseudomodes are always 
equivalent to discrete quasi modes, we can always obtain Markovian master equations for pseudomode-atom system. 



V. MARKOVIAN MASTER EQUATION FOR ATOM-DISCRETE QUASI MODES SYSTEM 



A key idea for treating the behaviour of a small system coupled to a structured reservoir is that although the behavior 
of the small system itself is non-Markovian, an enlarged system can obtained that exhibits Markovian dynamics — and 
which includes the small system, whose dynamics can be obtained later. In our example of a multilevel atomic system 
coupled to the the quantum EM field as a structured reservoir, we can proceed as follows. The overall system of the 
atom(s) plus quantum EM field is partitioned into a Markovian system consisting of the atom plus the discrete quasi 
modes and a flat reservoir consisting of the continuum quasi modes. The system Hamiltonian Hs is: 



Hs = ^ Vkhujk i^t^k ^ ^k^k) + ^ Tiuj^a\a^ 

k i 



k i 



whilst the reservoir Hamiltonian Hji is: 



Hr^ J dApc{A)hAb''{A)b{A) 



and the system-reservoir interaction Hamiltonian i?s-R is: 

Hs-R = Y. I ^Apc(A) \nW^(A)a^J^{A) + h.c. 



(95) 



(96) 



(97) 
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so that the total Hamiltonian is still equal to the sum of Ha, Hp and Haf, given in Eqs.(p^, (|l3|) and (|T5|). The 
distinction between the non-Markovian true mode treatment and the Markovian quasi mode approach is depicted in 
Fig. i 

It is of course the slowly varying nature of the coupling constants Wi{A) and the mode density Pc{A) which results in 
a Markovian master equation for the reduced density operator p of the atom-discrete quasi modes system. Rather than 
derive the master equation for the most general state of the reservoir, we will just consider the simplest case in which 
the reservoir of continuum quasi modes are all in the vacuum state. Again, the coupling constants Wi will be assumed 
constant so that no shift matrix Fij elements are present. The master equation is derived via standard proceedures 



(Born and Markoff approximations) |12 20|], which require the evaluation of two-time reservoir correlation functions 
in which the required reservoir operators are the quantities J dApc{A)Wi{A)b{A) and their Hermitian adjoints. To 
obtain Markovian behaviour we require the quantities pc{A)Wi{A)W* {A) to be slowly varying with A, so that the 
reservoir correlation time Tc (inversely proportional to the bandwidth of pc{A)Wi{A)W* {A)) is sufficiently short that 
the interaction picture density operator hardly changes during Tc- 

The standard procedure then yields the master equation in the Lindblad form: 



dp 
'dt 



Hs,P 



TTPcW^W* 



aj,pa. 



ajp, a\ 



(98) 



Direct couplings between the discrete quasi modes involving the Vij are included in the system Hamiltonian Us- 
Radiative processes take place via the atom-discrete quasi modes interaction also included in TJg, though still given as 
in Eq.(p^). The loss of radiative energy to the reservoir is described via the relaxation terms in the master equation. 
The diagonal terms where i = j describe the relaxation of the i th quasi mode in which the decay rate is proportional 
to Pc\Wi\'^. A typical decay rate Ti for the i th discrete quasi mode into the reservoir of continuum quasi modes will 
be: 



T,=2ttp,\W,\ 



(99) 



Note that the off-diagonal terms i ^ j involve pairs of discrete quasi mode operators dj and aj, so there is also a 
type of rotating wave approximation interaction taking place via the reservoir between these discrete quasi modes, as 
well as via direct Hamiltonian coupling involving the Vij . The standard criterion for the validity of the Born-Markoff 
master equation Eq. ( p8| ) is that Ttc <C 1 . Processes involving multiphoton excitation of the reservoir (such as may 
occur for excited multilevel atoms) can be studied using standard master equation methods, thereby enabling multiple 
excitation of the structured reservoir to be treated via the quasimode theory. 



VI. NON SLOWLY- VARYING MODE DENSITIES AND/OR COUPLING CONSTANTS 

The basic model treated in this paper is that of atomic systems coupled to a set of discrete quasi modes of the 
EM field, which are in turn coupled to a continuum set of quasi modes. Although expressions for the true mode 
coupling constant and the reservoir structure function have been obtained for the general case where the quasi mode 
density pc and the couphng constants Wi are not necessarily slowly varying functions of A (see Eqs.(^ and (pO|)), 
the usefulness of the results where this is not the case is somewhat limited. As indicated in the previous section, the 
master equation for the atom plus discrete quasi modes system will no longer be Markovian, so the enlargement of 
the system based on adding the discrete quasi modes to produce a Markovian system fails. 

Also, for the non slowly varying pc or Wi case, we can no longer link each discrete quasi mode to a pseudo mode. 
That this is the situation may be seen both from the general result for the reservoir structure function (Eq.(pO|)) or 
the specific result we have obtained for the case where there is a single discrete quasi mode (Eq.(p9|)). In the former 
case, the function Pn{'-^) would not be a polynomial of degree n, and therefore could have more than n roots, leading 
to more pseudo modes than discrete quasi modes. In the latter case involving just one discrete quasi mode, even 
having the mode density Pc{'^) (and hence r(cj)) represented by a single peaked function would result in D{lj) going 
from a single peaked function to a triple peaked function, corresponding to three pseudo modes. 

However, where pc or Wi are no longer slowly varying, an examination of the underlying causes for this variation 
may suggest replacing the present atom plus discrete and continuum quasi mode model by a more elaborate system 
that better represents the physics of the situation, but with only slowly varying parameters now involved. Fano 
diagonalisation based on such a more elaborate model could produce the desired link up with the pseudo mode 
approach and enable a suitable, enlarged system to be identified which has Markovian behaviour, as well as overcoming 
the problem of treating multiple reservoir excitations. One possible elaboration would be to add a further continuum 
of quasi modes that are fermionic rather than bosonic. 
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VII. CONCLUSIONS 



The theory presented above is mainly intended to apply to the important situation where the reservoir structure is 
actually due to the presence of a discrete system of quasi modes which arc coupled to other continuum quasi modes via 
slowly varying coupling constants. For example, the quantum EM field in high Q resonant cavities can be accurately 
described in terms of the quasi mode model which has these features, the discrete quasi modes being the cavity quasi 
modes (linked to the cavity resonances) with which the atoms inside the cavity interact, and the continuum quasi 
modes being the external modes. 

For this situation it has been shown that, for the present case of single quantum excitations, the pseudo mode method 
for treating atomic systems coupled to a structured reservoir of true quantum EM field modes can be obtained by 
applying the Fano diagonalisation method to the field described in an equivalent way as a set of discrete quasi modes 
together with a set of continuum quasi modes, whose mode density is assiimcd to be slowly varying. The interaction 
between the discrete and continuum quasi modes is treated in the rotating-wave approximation assuming slowly 
varying coupling constants, and the atomic system is assumed to be only coupled to the discrete quasi modes. The 
theory includes the true and continuum quasi mode densities explicitly. 

Expressions for the quasi mode operators Ui and 6(A) in terms of the true mode operators A{lv) (and vice versa) 
have been found, and explicit forms for the atom-true mode coupling constants have been obtained and related to 
the reservoir structure function that applies in pscudomodc theory. Wc have seen that the feature that leads to a 
pseudomode is the presence of a discrete quasi mode. Each discrete quasi mode corresponds to one of the pseudomodes, 
whose position Zi in the lower half complex plane is determined from the roots of a polynomial equation depending 
on the parameters for the quasi mode system. 

Although the behavior of the atom itself is non-Markovian, an enlarged system consisting of the atom plus the 
discrete quasi modes coiiplcd to a flat reservoir consisting of the continuum quasi modes exhibits Markovian dynamics, 
and the master equation for this enlarged system has been obtained. Using the quasimode theory, processes involving 
multiphoton excitation of the structured reservoir (such as may occur for excited multilevel atoms) can now be 
studied using standard master equation methods applied to the atom-discrete quasi modes system. Furthermore, 
cases with unusual forms of the reservoir structure function for single quantum excitation (for example, containing 
Lorentzians with negative weights) still result in Markovian master equations. Since for single quantum excitation 
the pseudomodes arc equivalent to discrete quasi modes, we can now always obtain Markovian master equations for 
pseudomode-atom systems via our approach. 

Although not so useful in such cases, the present theory does lead to general expressions for the true mode coupling 
constant and the reservoir structure function for single quantum excitation. These expressions are still valid for the 
general case whcirc; the; quasi mode density pc and the coupling constants Wi are no longer slowly varying functions 
of A. Howevcir. the mastcir equation for the atom plus discrete quasi modes system will no longer be Markovian, so 
the enlargement of the system based on adding the discrete quasi modes to produce a Markovian system fails. Also, 
for the non slowly varying pc or Wi case, we can no longer link each discrete quasi mode to a pseudo mode — there 
may be more pseudo modes than discrete quasi modes. In such cases it would be desirable to replace the present 
quasi mode system by a more elaborate quasi mode system involving only slowly varying quantities, and which better 
represents the underlying physical causes of the variation in Wj and Pc that occurs in the present model. This may 
make possible an extension of the Fano diagonalisation approach that still links quasi modes with pseudo modes, and 
results in a Markovian master equation for the enlarged atom plus quasi mode system. In such an elaborated system, 
the disadvantage of the present pseudo mode treatment in treating multiple excitations of the structured reservoir 
could still be removed. 

The treatment has been outlined in the case of a multilevel atom coupled to a structured reservoir of quantum EM 
field modes, but a similar approach would apply for any fermionic system coupled to a structured reservoir of bosonic 
oscillators. Extensions to fermionic reservoirs should be possible also. At present the treatment is restricted to cases 
where threshold and band gap effects are unimportant, but may be applicable to two-dimensional photonic band gap 
materials. Further extensions of the treatment to allow for atomic systems driven by single mode external laser fields 
are also possible, with the original atomic system being replaced by the dressed atom. 
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APPENDIX A: ATOMIC HAMILTONIAN 



As an example of writing the atomic Hamihonian in the form given in Eq. , consider a three level atom in a V 
configuration with upper states |2), |1) and lower state |0), whose energy is chosen for convenience to be zero. The 
atomic transition operators are ct^ = |2)(0| and = |1)(0| for the two optical transitions of frequencies UJ2 and wi, 
and a'^ = |2)(1| for the Zeeman transition of frequency 1^2 — wi- 

The form given in Eq.([lJ) is: 

Ha = rjihuJi {fJ^^i ~ ^i^t) + V2^i^2 {^2^2 ~ ^2 ^2) + V3^{^2 — wi) {^3^^ ~ '^3^3) 

= v^hu;^ - |0>(0|) + mni02 (|2)(2| - |0)(0|) + 773^^2 - ^i) (|2>(2| - |1>(1|) . (Al) 

This expression may also be written in the form 

Ha = nLUi\l){l\ + hLU2\2){2\ + hLu{\0){0\ + |1)(1| + |2)(2|), (A2) 

since by equating the coefficients of the three projection operators, we obtain a set of linear equations for the ?7i, 772, Tj^ 
and LO which are solvable — in fact the solutions are not even unique. These equations are: 

772^^2 + m{^2 - C^l) = 0^2 + W 
IJiLJi — rj'i{uJ2 — OJl) = UJi+W 

-rjiUJi - 772072 = 0;. (A3) 
Adding these equations and then substituting into the first two gives: 

w = -i(a;i+W2) (A4) 

+773(^2 - uji) ^ -^^1 + ~ ^2)^2 
2 1 

-?73(w2 - ^i) = - - 2'^2. (A5) 

The last two equations do not produce a unique solution for 7/1,772,773. We can arbitarily choose r/3 — for the low 
frequency transition, and then we find that: 

1 1 (CJ2 - ^1) 

3 6 oji 

V2 = \ + \^-^^^^. (A6) 

O O UJ2 

This gives »7i = 772 = for two degenerate optical frequency transitions. 

Comparing the two expressions f or H a in Eqs.(Al), (|A^), where 771,772 are given by Eq.(|A(]|) (with 773 set to zero) 
and u! by Eq.(A4), we see that Eq.(A2) gives the atomic energy apart from the constant term — ■i?7(a;i + CJ2). 



APPENDIX B: SCALING FOR MODE DENSITIES p(a;),pc(A) EQUAL TO UNITY 



The equations presented in the first part of Section [II are based on true and quasi continuum mode densities that 
are not necessarily equal to unity. To compare our expressions with those in the Ref. fl^ , we now set out the scalings 
needed for the various quantities to give the Hamiltonians equivalent to Hp and Haf in either true or quasi mode 
forms (Eqs.(^6|), (|l3|), ( |l8|) and ([Tsl)) in which the mode densities p and pc are made equal to unity. The creation 
and anihilation operators are no longer dimensionless, the coupling constants and angular frequencies do not have 
dimensions of frequency and the expansion coefficients are not dimensionless. The scaled quantities appearing in the 
Hamiltonians or relationships between annihilation operators will be denoted with a superscript '^^K 

The following replacements made to the annihilation and creation operators: 



v/^6(A)^fe(^)(A) (Bl) 

^^:(A)6t(A) ^ fe(^)t(A) (B2) 

^M^A{iu) ^ A^'\iu) (B3) 

V^it(^)_i(^)t(^)^ (B4) 
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to the coupling constants: 

./^)W,{^)^wl'\A) (B5) 

and to the expansion coefficients : 

V^a.(a;) ^ a,{ujf'^ (B6) 
v/pe(A)pH/?(c.,A)^/3(^)(c^,A) (B7) 

wiU give the Hamiltonians equivalent to Hp and Haf in either true or quasi mode forms (Eqs . ([l6|) , , ([l8|) and 
(p^)) in which the mode densities are put equal to one. In addition the modified forms of the relationships between 
true and quasi mode annihilation operators (Eqs.(p2|) and (^2[)) can be obtained in which p and Pc are made equal to 
unity, as can the revised forms of the commutation rules. The latter are: 

[^(-)(A),&(«)t(A')] =5(A_ A') (B8) 
[A^'\uj),A^'^\uj')]^5{uj~Lo'). (B9) 

In addition the various equations for the Fij{uj),ai{Lo), (3{uj, A),z{uj), g'' {uj) , D'^ (uj) now apply with p and pc put equal 
to unity. It should be noted that the quantities ai,uji,Vij,A,u} are not replaced, nor are any of the atomic quantities 



APPENDIX C: THE HAMILTONIAN Hp IN DIAGONALISED FORM 

We show by starting with the field Hamiltonian in the quasi mode form Eq. (p^) , substituting the solutions for 
ai{u}) and /3(a;, A) into the expressions for di and &(A) given in Eqs.jS^) and then evaluating the result, that the 
Hamiltonian in the true mode form Eq.(|6|) is obtained. The symmetry conditions Fij — F*^ and Vij — v*^ are used 
throughout. 

Using the expressions for hi and 6(A) given in Eqs.(^2|) the Hamiltonian in the quasi mode form Eq.(^3|) is then 
given by: 



Hp = h j dujp{uo) j dJ p{J)A\uj)A{J)I{uj,J) (CI) 

where the function is: 

/(w,w/) =^a;,a,(a;)a*(w/)+ f dAp,{A)Af3{A,Lu)f3* {A,lu^) 

+ J2 ''^3(^^i^h^i^^)+Y. J dAp,{A)W,{A)a,iu)P*{A,J) 

+ rfAp,(A)W,*(A)ar(a;/)/3(A,a;) (C2) 



Substituting for /3(w, A) in terms of the ai{uj) from Eq.(|25|), using the expression ( pTj) for Fij and then Eq. (|26|) for 
the ai{u!) we get for certain contributions within the last two terms in Eq.(p2|): 



dAp,iA)W.iA)P*iA,u;/)=-icj,~u;/)a*icj/)~ ^ J^*.a*iu;/) 
j dAp,{A)w:{A)l3{A,uj) =-{uj,-Lo)a,{Lo)- ^ J^,»a,(w) (C3) 

leading to 

^ j dAp,{A)W,{A)a,{uj)P*{^,^') = - - - E (C4) 

E j dAp,{A)W:{A)a*{J)P{A,u) = - ^(c., - u;)a.(c.)a*(^/) - ^ v,,a,{u)a*{J) (C5) 
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In the second term of Eq.(|C2|) substitution for l3{uj, A) and f3*{ui/ ^ A) in terms of the ai{u}) and a*(w/)from Eq.(|2E 
and then using the Eqs.(|37|) for manipulating principal integrals and delta functions leads to: 

y dAp,(A)A/3(A,c^)r (A,c^/) = ^{[y dAp,(A) x 



xA 



W,{^)a,{Lo)Wn^)a*Auj')] 



uj/ — uj \uj — A cj/ — A ' i\ J i.\ J ] 
TT^5{ijj - uj')pc{uj)ojW,{uj)a,{uj)W*{ijj)a*{uj) 
V 



UJ — uj/ 
V 



2 \^ J^] < 

p,{J)z*{uj/)W,{uj/)a,{uj)W*{J)a*AJ) 



uj/ — UJ' 



p,{uj)z{uj)W^{uj)a,{uj)Wnuj)a*{uj') 



3 V~ /"J > 

/) 



+ujS{uj — uj/)pciuj)z{uj)z*{uj)Wi{uj)ai{uj)W*{uj)a*{u>)} 



Then using Eq.(p8|) we show that 

A 



V I V 

UJ' 



UJ - A uj/ - A J \ UJ - A uj/ - A 
and following the introduction of the Fij from Eq.(|27|) we get: 

/ dAp,(A)A/3(A,c^)/3*(A,c^/) = y^{uj^^Fj,{uj)aMa*{J) 

J ^—^ UJ' — UJ ■' 



(C6) 



(C7) 



^UJ- 



Uj/ ~ UJ 

V 



F,,{J)a,{Lo)a*(uj/) 



I _ p,[uj)z{uj)W,{uj)a,{uj)W*(uj)a*(LO^) 



UJ' — UJ 

I V 



■p,{J)z*{J)W^{J)a^{uj)W*{J)oL*{uj/) 



+ujb{uj - J)p,{lo){t^ + |z(w)|2]Ty,(w)a,(w)W*(cu)a*(cj)} 



(C8) 



The last term is just ujb(uj — uj/)/ p(uj) using the normalization condition Eq.(|39|), the Pc(^) factor cancelling out.. 
The next step is to eliminate the Fy using Eq.(p6|) for the ai{uj) twice. After further algebra using Eq.(^) again we 
find that: 



dAp,(A)A/3(A, w)/3*(A,u;/) = uj5{uj - uj/)/p{uj) + ^ 

i 

-{uJ+J)Y.a,{uj)a*{J)+ ^ v*M'^)a*{uj/) 



(C9) 



The results in Eqs.(C9), (C4) and (C5) can be substituted back into Eq.(C2) for I{uj,uj/). It is found that there is 
extensive cancellation leading to the final expression: 



I(uj, uj/) = uj5{uj — uj/)/p{uj) 



and hence the Hamiltonian Hp in Eq.(Cl) is now in its true mode form: 

Hf ^ dujp{uj)fiujA''{uj)A{uj) 



(CIO) 



(CU) 



thus showing that the true and quasi mode forms of Hp are equal. 
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Reservoir 

Illustration of a two- level atom coupled to a structured reservoir. 
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FIG. 2. Three pictures of the coupled atomic system. In the true-mode picture the atom is coupled directly to true modes 
which have structure. In the quasi-mode picture the atoms are coupled to quasi-modes which are in turn coupled to external 
quasi-modes. In the pseudo-mode picture the atoms are coupled to dissipative pseudomodes. 
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